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Abstract 

Anisotropic mesh adaptation has been successfully applied to the numerical solution of 
partial differential equations but little considered for variational problems. In this paper, 
we investigate the use of a global hierarchical basis error estimator for the development of 
an anisotropic metric tensor needed for the adaptive finite element solution of variational 
problems. The new metric tensor is completely a posteriori and based on residual, edge 
jumps and the hierarchical basis error estimator. Numerical results show that it performs 
comparable with existing metric tensors based on Hessian recovery. A few sweeps of the 
symmetric GauB-Seidel iteration for solving the global error problem prove sufficient to 
provide directional information necessary for successful mesh adaptation. 
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1. Introduction 

Variational problems arise from many areas of science and engineering and have a nat- 
ural formulation with which the governing equation can be derived through minimization. 
Often, a variational problem can be transformed into a boundary value problem of partial 
differential equations (PDEs) and solved by methods specially designed for PDEs. Unfortu- 
nately, these methods are not designed to take structural advantage of variational problems 
and many researches argue that the variational formulation should be used as a natural 
optimality criterion for mesh adaptation; e.g., see [3llll[5]. In this article, we consider mesh 
adaptation based on the variational formulation. 

Most of the existing work for the adaptive numerical solution of variational problems 
employs isotropic mesh adaptation where the size of mesh elements varies from place to 
place according to a certain local error estimate while their shape is kept close to equilateral. 
On the other hand, the ability to allow the size, shape and orientation of mesh elements to 
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vary can significantly improve the accuracy of tlie solution and enliance tlie computational 
efficiency; this is especially true for problems exhibiting distinct anisotropic features. While 
it has attracted considerable attention from many researchers and been successfully applied 
to the numerical solution of PDEs, anisotropic mesh adaptation has rarely been employed 
for variational problems, especially when combined with a posteriori error estimates. 

The objective of this article is to investigate the application of a hierarchical basis a pos- 
teriori error estimator (HBEE) to anisotropic mesh adaptation for variational problems. We 
use the M -uniform mesh approach pL5j| where an adaptive mesh is generated as a uniform 
mesh in the metric specified by a symmetric and strictly positive definite tensor, M. The 
key to the success of the approach is to define a proper metric tensor. Recently, Huang and 
Li |TH] developed a metric tensor for the adaptive finite element solution of variational prob- 
lems. In the anisotropic case, M is semi-a posteriori: it involves residual and edge jumps, 
both dependent on the computed solution, and the Hessian of the exact solution, which 
provides necessary directional information for anisotropic mesh adaptation and is typically 
approximated using a Hessian recovery technique [18]. On the other hand, an HBEE was 
used by Huang et al. [IT] to develop a metric tensor for the adaptive finite element solution 
of the boundary value problem of a second-order elliptic differential equation. In this paper 
we combine the approaches in [TTJ [18] and develop a metric tensor for variational problems 
based on an HBEE and the underlying variational formulation. The obtained metric tensor 
is a posteriori in the sense that it is based solely on residual, edge jumps, and an a posteriori 
error estimate. This is in contrast to most previous work where M depends on the Hessian of 
the exact solution and is semi-a posteriori or completely a priori; e.g., see [3 |H1 El HH IISl [IE] • 
Numerical results show that the new metric tensor works well for all numerical examples 
and is comparable in performance with those based on Hessian recovery. 

The outline of the paper is as follows. Section [2] is devoted to the description of a 
general variational problem, its finite element approximation, and the general procedure of 
the M-uniform mesh approach for mesh adaptation. In Section [3} a bound is derived for 
the variation of the underlying functional. It involves residual, edge jumps, and an HBEE 
and is minimized on M-uniform meshes for an optimal M. Numerical examples are given 
in Section |4j Finally, concluding remarks are given in Section [5] 

2. Finite element approximation and mesh adaptation 

Consider a general functional of the form 



where -, ■) is a given smooth function, f2 C M'^ (c? = 1, 2, 3) is the physical domain and 
Uq is a properly selected set of functions satisfying the homogeneous Dirichlet boundary 
condition 
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v{x) =0 Va; G dVL. 
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In this paper, we consider only the homogeneous Dirichlet boundary condition. Other 
boundary conditions can be either converted to the homogeneous Dirichlet boundary con- 
dition [I3l or treated without major modifications. The corresponding variational problem 
is to find a minimizer u ^ Uq such that 

I[u] = mm I [v]. (2) 

We assume that F satisfies suitable conditions so that the variational problem ^ has a 
unique solution. We also assume that a linearization of the first variation of the functional 
(i.e. B[u; ■, ■], see Q below) is convex around the solution of the minimization problem. 
The latter is needed for the error problem ^ defined in Section 3.2 to have a unique 



solution and for its GauB-Seidel solution to converge. However, we do not assume that F is 
quadratic although some theoretical analysis in this paper is valid only for convex quadratic 
functionals. 

A necessary condition for u to be a minimizer is that the first variation of the functional 
vanishes. This leads to the Galerkin formulation 

6I[u,v]= / {Fu{x,u,'Vu) V + F^u{x,u,'Vu) ■ Vv) dx = 0, MveUq 
Jn 

where Fu and Fy^ are the partial derivatives of F with respect to u and Vm, respectively. 
The linearization of the first variation 6I[-, ■] around the minimizer u with respect to the 
first argument is given by 

f f d'^F d^F d^F 

B[u; w,v]= v—^{x, u, Vu)w + v {x, u, Vu) ■ Vw + w {x, u, Vm) ■ Vf 

d^F \ 
+Vw-———{x,u,Vu)-Vv]dx, ^w,veUo. (3) 



dVudVu 

Given a triangulation Th for fi, let Uq C Uq be the associated linear finite element space. 
Then, a linear finite element approximation of u is defined as Uh G Uq such that 



I[uh] = min I[vh\. 

The solution Uh can also be found by solving the corresponding Galerkin formulation 

5I[uh, Vh] = / {F^{x, Uh, Vuh) Vh + Fv«(a;, Uh, Vuh) ■ Vvh) dx = 0, \/vh G Uq. 
Jn 



(4) 



For the adaptive finite element solution, the mesh Th is generated according to the behaviour 
of the error of the approximation Uh- We follow the so-called M-uniform mesh approach 
[15] with which an adaptive mesh is generated as a uniform mesh in the metric specified by 
a symmetric and strictly positive definite tensor M = M{x). Such a mesh is called an M- 
uniform mesh. A scalar metric tensor will lead to an isotropic mesh while a full metric tensor 
will generally result in an anisotropic mesh. In this sense, the mesh generation procedure is 
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the same for both isotropic and anisotropic mesh generation. The key to the approach is to 
choose a proper metric tensor M. 

Once a metric tensor M has been chosen, a sequence of mesh and corresponding finite 
element approximation are generated in an iterative fashion. We start with an initial mesh 
rl:'\ For every mesh T^*-*'' we solve the variational problem with (j) for finite element 

approximation uf* . The new mesh T^*'*^^^ is generated as an almost M-uniform mesh in the 
metric specified by M^*^ which is computed based on u'^ and 7^^'\ This yields the sequence 

/^-(O) rrh N (0) , .(0) .^(1) jjh ^ (1) ,^(1) 

The process is repeated until a good adaptation is achieved. 

In the following section we derive the metric tensor M for use in anisotropic mesh adap- 
tation for the variational problem (|2|. First, we derive a bound on the variation 5I[uh,v] 
and then define M such that the obtained bound is minimized on corresponding M-uniform 
meshes. 



3. Monitor functions based on hierarchical basis error estimation 

3.1. Error bounds for variational problems 

Let eh = u — Uh he the error of the finite element solution Uh- It is known [12j that for 
a uniformly convex quadratic functional in the form ([T]), there exists a positive constant /3 
such that 

||Ve,||i.(^) </3|5/[e;„e,]| =/3|5/[n,,e,]|. 
In this case, the quantities 

l^/K,e.]| and f^I^^V (5) 

are equivalent to || Ve/i||^2(-Q), and minimizing their bounds is equivalent to minimizing error 
bounds. Consequently, it is reasonable to define M based on bounds for these quantities 
even for other functionals. 

We begin with deriving a bound for e/j] |. For simplicity we denote 

Fu{x) = Fu{x, Uh{x), Vuh{x)), F^u{x) = F^^{x, Uh{x),Vuh{x)). 

Then for any v eUq 



6I[uh,v\ 



/ [Fu{x, Uh, Vuh) V + F^u{x, Uh, Vuh) ■ Vv] dx 

V / [F^{x)v + Fy^{x)-Vv]dx. (6) 



From the divergence theorem, (ml) can be rewritten as 



6I[uh,v]=J2 / [Fuix) - V ■ Fy^{x)] V dx + J2 / vF^uix)-nds, (7) 



where n denotes the outward unit normal to the face dK. 

Now, let dTh be the collection of all faces of mesh Th and K and K' be the two elements 
sharing a common face 7. We define residual r/j and edge jump Rfi as 



(F^uix) ■ n^)\K + (F^uix) ■ n^)\K' a? G 7 ^'j E dTh\dn, 
a; G 7 V7 G 



Rh(x) 



Then becomes 

6I[uh,v]= Y / rh{x)vdx+ ^ / Rh{x)vds. 

Taking v = = u — Uh va. the above equation and using Schwarz's inequality, we obtain 

\5I[uh,eh]\< Yj / \rh{x)eh\dx+ ^ / \Rh{x) eh\ds 
KeTh -yean •^'^ 

- 'Y W^f^^L^iK) \\'^h\\L^(^K) + 11-^/11^2(7) lk/i|lL2(7) 



(8) 



KeTh 



\^h\\L2^K) \\eh\\L2^K) + - ^ I|-R/i|Il2(7) II'^^IIl2{7) 

■yedK 



For further derivation, we need to estimate ||eh||2,2(^) and ||e/i||j;^2(^). In [i8\, the authors 
employ anisotropic interpolation error bounds on elements and element faces. This results 
in a semi-a posteriori metric tensor which involves residual and edge jumps, both dependent 
on the computed solution, and the Hessian of the exact solution. In the following, we use 
an a posteriori hierarchical basis error estimate. 

3.2. An a posteriori error estimator based on hierarchical bases 

The computation of the error estimator is based on a general framework, details on which 
can be found among others in the work of Bank and Smith |2] or Deufihard et al. [10] . The 
approach is briefly explained as follows. 

Recall that Uh G f/g is a linear finite element solution of the Galerkin formulation Q and 
the error is Ch = u — Uh. Consider the extension of Uq to a subspace of piecewise quadratic 
functions with the linear span of the edge bubble functions W'^, i.e. 

Uo ul; © w^. 
5 



Recall also that 



6I[uh + eh,v] =0 Vt; G Uq. 
The error estimate Zh is then defined as the solution of the approximate linear error problem 

{Find Zh G W'^ such that 
6I[uh, Wh] + B[uh; Zh, Wh] = ^Wh e 

Recall that B[uh',-,-] is the bilinear form resulting from the linearization of ^/[■, ■] about 
Uh with respect to the first argument. The estimate Zh can be viewed as a projection of 
the error onto the subspace W'^. By construction, HhZh = for the linear finite element 
interpolation operator 11^. Thus, Zh — ^h^h = Zh- 

This definition of the error estimate is global and its solution can be costly. To avoid 
the expensive exact solution in numerical computation, we employ only a few sweeps of the 
symmetric GauB-Seidel iteration for the resulting linear system, which proves to be sufficient 
for the purpose of mesh adaptation (see [17] or numerical examples in Section 111). 

It is noted that under the convex assumption of B[u]-,-], the error problem (|9| is solvable 
and has a unique solution when Uh is sufficiently close to u. Moreover, Gaufi-Seidel iteration 
is convergent when applied to the solution of 

3.3. Optimal metric 

We now use the error estimate Zh to develop the metric tensor. Instead of directly using 
the bound (Isl), we use 



E[uh,Zh] = 

K&Th 



\''^h\\L2(K) \\^h\\L2(K) + 2 ll-^/illL2{7) lk/illL2(^) 



(10) 



which is obtained by replacing Ch in ([8j) with Zh- Theoretically, it is unclear if the bound (|8 
is bounded by E[uh,Zh] although our numerical results show that the metric tensor based 
on E[uh, Zh] leads to proper mesh adaptation for all examples we have tested. On the other 
hand, this is somewhat related to the saturation assumption, which basically states that 
the quadratic finite element approximation provides a truly better approximation than the 
linear one. The saturation assumption is known to be valid for some cases under suitable 
conditions for isotropic meshes [H]. It is still unknown if such results hold for anisotropic 
meshes. 

Recalling that HhZh = and using element-wise interpolation error estimates in p^, we 
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have 



lk/i|lL2(A') — W^h — ^hZh\\L2(^x-j 



Ik 



(tr 



E l II ||2 
I I II^/iIIl2(7) 



< C 

= C\K\hT{{F'^f\HK 
1 

M 



{{F-f\H{z,)\F'^)Ydx 



(11) 



< C 



^ 1 



element i^, and C is a 



where tr(-) denote the trace of a matrix, Hxizh) is the Hessian of on i^', = 
^/Wj^izh), \K\ is the volume of K, F^ is a mapping from the reference element K to 
element i^, and C is a constant independent of Th and Zh. Thus 

^ L2(7)= 5Z l7lMl^h|lL2(7)I 



■y£dK 



' '""11^2(7) 
1 



5Z U| 11^^11^2(7) 



V7ee-ft: 

Substituting (11) and (12) into (10) leads to 



<c( E I7NI^/.IIl^(7)) tr((F;,f 



(12) 



KeTh 
KeTh 



1^1^ lk/^llL^(x) + E 1^1' II^'^IIl2(7) ((f;,)^ f;,) 

yedK 

I^^i^i^l |Ait.((F;,)^|ff*-(^-,.)|Fi) 



where (■) denotes the scaled norm 



' ' -yedK 



- I ,1/2 II'"IIl2(u;) - 



We now use bound (13) to define the metric tensor M. To ensure that M is strictly 
positive definite, we first regularize the bound with a positive constant au (to be determined), 
i.e., 



K&Th 



X 



Ken 



X 



ah + {rh)L2(K) + 5Z 1^1 
\K\tT{{F'j,f{aHl+\HK{zh)\)F'j,) 

K\tT{iF'j,fHKA^h)F'K) 



L2(7) 



L2(7) 



■y&dK 



(14) 



where 



ah 



The optimal metric tensor is obtained by minimizing bound (|14|) for M-uniform meshes. 

(15) 



It is known [16] that an M-uniform mesh satisfies the alignment condition 

1 



d 



tr {{F'^fMKF'^) = det {{F'^YUkF'^) 



and the equidistrihution condition 



Pk \K\ 



(16) 



where = ^^eTh Pk\^\^ is an average of M over K, px = ■\/det(M/f), and the 
number of elements of Th- 

We now pay our attention to the tr(-) factor in (14). Notice that in general, 

1 



d 



tr {{F'^fHKAzh)F'K) > det {{F'j,fHKA^h)F'Ky. 



From (15) we can see that the equality in the above inequalities holds if we choose M = 
in the form 

Mk = OkHkA^u) VK e Th (17) 

for some scalar function 9 = Ok- Indeed, with this choice of Mk the alignment condition 
(15) reads as 



^tr {{F'j,fHK,aF'K) = det {{F'j,f HK,aF'Ky = det {HKA^h)^ 



(18) 
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where we have used |det(F^)| = \K\ and assumed \K\ = 1. Substituting (18) into (14) 
yields 



E[uh,Zh]<Cal J2 



'L2(7) 



d+2 



X \K\ ^ det {HK,a{zh)y ■ 



(19) 



Next, we use the equidistribution condition ([16j) to determine 9 = 9^ in (17). From Holder's 
inequahty, we have 



d 
d+2 



or 



K \ K 

E d+2 d+2 2 

i\K\pK) ' > K) ' 

K 



Oji 

N 



(20) 



with the lower bound being attained for a mesh satisfying (16). Comparing the left-hand 
side of (20) with the right-hand side of (19) suggests that p = px he defined as 

det {HK,aiZh))^' 

n. I ' — ' 

d 

-y 1 d+2 



Pk 



d 
d+2 



1 



det [I + — \HK{zh) 

Olh 



d+2 



(21) 



From relations px = \/ det(M;^) and Mk = 9KHK,a{zh) we can obtain 9k- The metric 
tensor Mk is then given by 



(22) 









(l+ — \HKizH)\' 


V 




\ ah J 




Olh 



With this choice of pk (and Mk)-, the right-hand side of (19) attains its lower bound for 
a mesh satisfying the equidistribution condition (16). Then, the variation of the functional 
([T| has an upper bound as 



d+2 2 



E\uh,Zh\ < Cai{ah)^N—^. 

To complete the definition of the metric tensor, we need to choose the regularity param- 
eter ah- Following [15], we choose it such that 

crh=J2 Pk\^\ = 
KeTh 

9 



or 



Ken 



(Rh) 



L2(7) 



d 
d+2 



det 



J + — \HK{Zh) 



d+2 



2\VL\ 



(23) 



With this choice, roughly fifty percents of the mesh elements will be concentrated in the 
regions of large p [12]. It is easy to show that (23) has a unique solution since its left-hand 



side is monotonically decreasing with ah increasing and tends to |n| as — )■ oo and to -|-cxo 
as ah — 0. Moreover, it can be solved using a simple iteration scheme such as the bisection 
method. 

Once ah is computed, the adaptation function p and the metric tensor M can be de- 



termined by (21) and (22), respectively. A new mesh can then be generated based on the 
metric tensor. 

This definition of the metric tensor involves the residual r/i, the edge jump i?^, and the 
HBEE Zh- AH these quantities are based on the computed solution; in this sense, (22) is 
a posteriori. 



4. Numerical results 

In this section, we present some numerical results for a selection of two-dimensional 
problems with an anisotropic behaviour and compare the following metric tensors: 

• isotropic a posteriori metric tensor [18] based on the variational formulation, residual, 
and edge jumps. 



• anisotropic semi-a posteriori metric tensor based on the variational formulation, 
residual, edge jumps, and Hessian recovery, 

• anisotropic a posteriori metric tensor based on the variational formulation, residual, 
edge jumps, and hierarchical basis error estimator (HBEE) developed in the previous 
section. 

The HBEE is computed as a global error problem. A few symmetric Gaufi-Seidel itera- 
tions are employed to obtain an approximation for Zh until the relative difference of the old 
and the new approximations is under a given tolerance GS-RTOL = 0.01. More details on 
several implementation issues such as mesh quality measure and computation of the error 
estimator can be found in [TTJ Section 4]. 

4.1. A quadratic functional 

As a first example we consider a quadratic functional 



I[u] 



^ \Vuf -uf 
10 



dx (24) 




(a) Example 4.1 the exact solution. 



iso + variational (residuals) 
aniso + variational (residuals + HBj 
aniso + variational (residuals + Hess, recovery) 
aniso (HB error estimator) 




number of elements 



(b) Example 



4.1 



error comparison. 



4.1 



Figure 1: Example 
II Ve/i||^2(t^) against the number of elements 



surface plot of the exact solution, (b) the finite element error 



with Q = (0, 1) X (0, 1). The function / and the Dirichlet boundary conditions are chosen 
such that the exact solution of the minimization problem ^ is given by 

Uexact{x, y) = tanh(60x) — tanh (60(x — y) — 30) . 



A surface plot of Uexact is given in Fig. [Ta) The solution exhibits a strong anisotropic 
behaviour and describes the interaction between a boundary layer along the x-axis and a 
shock wave along the line y = x — 0.5. 

Note that this variational problem is equivalent to the boundary value problem 

-Au = f in Q, 

U = Uexact OU 8^1. 

Thus, we also compare metric tensors based on the variational formulation described above 
with the one based on the hierarchical basis error estimator and PDF developed in [T7] . 

shows II Ve/i||j;^2(Q) (the finite element solution error measured in the semi- 



Figure 



lb 



norm) agamst the number of elements for all metric tensors; adaptive mesh examples as well 
as close-up views near the point (0.5,0), where the boundary layer meets the shock wave, 
are given in Figs. [2]and[3j 

First, we observe that all methods provide good mesh density, but anisotropic adaptation 
(Fig. |3| provides much better mesh alignment with the boundary layer and the shock wave 
front than isotropic adaptation (Fig. |2]). Once again, it shows that anisotropic adaptation 
is advantageous and should be favored for such problems. Meshes and the solution errors 
obtained with anisotropic adaptation are very similar. Moreover, mesh adaptation by means 
of the variational formulation and HBEE is comparable with those obtained with both the 
Hessian recovery-based method in [18j and the PDF HBEE method in fTTJ. 
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Figure 2: Example 4.1 (isotropic): adaptive mesh and a close-up view at (0.5,0) obtained 
with the isotropic metric tensor based on variational formulation with residual and edge 
jumps; 664 vertices and 1226 triangles, ||Ve/i||^2 = 2.8, maximum aspect ratio 3. 



One may observe that the PDE HBEE method produces a slightly smaller error. The 
reason for this behaviour can be attributed to the use of the residual and the edge jump in 
the metric tensor based on the variational formulation: it can be seen as a data smoothing, 
which prevents the sharp change of the element aspect ratio. The effect can be observed in 
Fig. [3] the maximum aspect ratio for the metric tensor based on the residual and HBEE 
is 33 (Fig. 3b), whereas the maximum aspect ratio of the mesh obtained by means of just 
HBEE is 45 (Fig. 3c). Hence, the PDE HBEE based metric tensor might be a better choice 



for quadratic functionals with strong anisotropic features. 

4-2. A non- quadratic functional 

The second example is an anisotropic variational problem [6] defined by the functional 



I[u] 



(1 + |Vm|^)^^ViOOOmJ 



dx 



(25) 



with Q = (0, 1) X (0, 1) and the boundary condition 



u 



u 



1 on X = or X = 1, 

2 on?/ = Oor?/ = l. 



Unlike the previous example, the functional (25) is not quadratic. Hence, the quantities 

\5I[uh,eh] \ and (\5I[uh,eh] \ / \ 

^^'ilL2(r2) ) mathematically equivalent to || Ve/i||^2(Q); 

this example is a good test for the formulas of the metric tensor based on the variational 
formulation. Analytical solution is not available for this example. A computed solution in 
Fig. |4] shows that the mesh adaptation challenge for this example lies in the resolution of 
the sharp boundary layers near x = and x = 1. 

Adaptive meshes obtained with the three different metric tensors (isotropic, Hessian 
recovery-based, and HBEE-based) are given in Fig. Isl They all have correct mesh concen- 



tration, but the anisotropic metric tensors (Figs. 5b and 5c ) provide a much better alignment 
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0.5 1 0.5 0.6 0.7 

(a) Variational formulation with residual, edge jumps, and Hessian recov- 
ery: 663 vertices and 1235 triangles, |lVe;i|jj^2 = 6.2 x 10~^, maximum 
aspect ratio 34. 





(c) HBEE: 677 vertices and 1256 triangles, |lVe,i||^2 = 6.0 x 10"\ maxi- 
mum aspect ratio 45. 



Figure 3: Example 4.1 (anisotropic): adaptive meshes obtained with different anisotropic 
metric tensors (left); close-up views at (0.5,0) (right). 
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Figure 4: Example 



4.2 



surface plots of the numerical solutions. 



with the boundary layers. Again, both anisotropic meshes are comparable, although mesh 
elements near the boundary layer in the HBEE-based adaptive mesh have a larger aspect 
ratio than elements of the mesh obtained by means of the Hessian recovery. This could be 
due to the smoothing nature of the Hessian recovery: usually, it operates on a larger patch, 
thus introducing an additional smoothing effect, which affects the grading of the elements' 
size and orientation. The global hierarchical basis error estimator does not have this hand- 
icap and, in this example, the mesh obtained by means of HBEE is slightly better aligned 
with the steep boundary layers. 

4-3. Image processing 

Our third example is an energy functional 

I[u]= [ \{p-uf + {l + \VufY^^]dx, (26) 

which is used in image processing with the observed image p and the reconstructed image u 
[T]. In our computation, we choose Q = (0, 1) x (0, 1), 

p = 1/ (l + gl000(a'+y-1.25)^ 

and the boundary condition 

u = p on dVt. 

The analytical solution is not known; a computed numerical solution is shown in Fig. [6] and 
Fig. [7] shows examples of adaptive meshes. 

As in the previous example, the anisotropic metric tensors are comparable and provide a 
better mesh adaptation than the isotropic one. Again, the HBEE-based mesh has a slightly 
larger maximum aspect ratio. 
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(a) Variational formulation with residual and edge jumps (isotropic): 644 
vertices and 1097 triangles, maximum aspect ratio 3.4. 





0.5 1 0.05 0.1 

(b) Variational formulation with residual, edge jumps, and Hessian recov- 
ery (anisotropic): 649 vertices and 1160 triangles, maximum aspect ratio 
15. 





(c) Variational formulation with residual, edge jumps, and HBEE: 639 
vertices and 1143 triangles, maximum aspect ratio 51. 



Figure 5: Example 4^ adaptive meshes obtained by means of isotropic, anisotropic Hes- 
sian recovery-based, and anisotropic HBEE-based metric tensors (left); close-up views at 
(0,0) (right). 
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Figure 6: Example 



4.3 



surface plots of the numerical solution. 



5. Concluding remarks 

Numerical results confirm the conclusion of [TTj that a global HBEE can be a successful 
alternative to Hessian recovery in mesh adaptation; a fast approximate solution of the global 
error problem is sufficient to provide directional information for anisotropic mesh adaptation. 
They also confirm the conjecture that good mesh adaptation does not require a convergent 
Hessian recovery or an accurate error estimator, but rather some additional information of 
global nature, although it is still unclear which information exactly is necessary for successful 
anisotropic adaptation. 

For quadratic functionals, numerical results suggest that algorithms which involve non- 
directional information such as residual sometime can be slightly disadvantageous for prob- 
lems with sharp anisotropic features. Such algorithms can smooth out the metric tensor 
and therefore prevent the rapid change of element size and orientation, which is highly de- 
sired for this class of problems. Also, as mentioned in [17], Hessian recovery could result in 
unnecessarily high mesh density for problems with discontinuities. Global hierarchical error 
estimation seems less affected by these issues and, depending on the underlying problem, 
can provide a more robust solution. 
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0.5 1 ' 0.9 0.95 1 

(a) Variational formulation with residual and edge jumps (isotropic): 657 
vertices and 1207 triangles, maximum aspect ratio 2.7. 




0.5 1 0.9 0.95 1 

(b) Variational formulation with residual, edge jumps, and Hessian recov- 
ery (anisotropic): 662 vertices and 1219 triangles, maximum aspect ratio 



3.5. 





0.5 1 0.9 0.95 1 

(c) Variational formulation with residual, edge jumps, and HBEE: 656 
vertices and 1209 triangles, maximum aspect ratio 4.7. 



Figure 7: Example |4.3[ adaptive meshes obtained by means of isotropic, anisotropic Hes- 
sian recovery-based, and anisotropic HBEE-based metric tensors (left); close-up views at 
(1,0.25) (right). 
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